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Abstract 

This is a proceedings article reviewing a recent combinatorial construction of the su(n)fe 
WZNW fusion ring by C. Stroppel and the author. It contains one novel aspect: the exphcit 
derivation of an algorithm for the computation of fusion coefficients different from the Kac- 
Walton formula. The discussion is presented from the point of view of a vertex model in 
statistical mechanics whose partition function generates the fusion coefficients. The statistical 
model can be shown to be integrable by linking its transfer matrix to a particular solution of the 
Yang-Baxter equation. This transfer matrix can be identified with the generating function of 
an (infinite) set of polynomials in a noncommutative alphabet: the generators of the local affine 
plactic algebra. The latter is a generalisation of the plactic algebra occurring in the context 
of the Robinson-Schensted correspondence. One can define analogues of Schur polynomials in 
this noncommutative alphabet which become identical to the fusion matrices when represented 
as endomorphisms over the state space of the integrable model. Crucial is the construction of 
an eigenbasis, the Bethe vectors, which are the idempotents of the fusion algebra. 



§ 1. Introduction 



Wess-Zumino-Novikov-Witten (WZNW) models are an important class of confor- 
mal field theories (CFT) distinguished by their Lie algebraic symmetry. Due to this 
symmetry the primary fields of WZNW theories are in one-to-one correspondence with 
the integrable highest weight representations of an affine Lie algebra; see e.g. the text 
book [Ij for details and references. Consider the su{n)k WZNW model, then the set of 
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all dominant integral weights of level k G Z>o is given by 

n 

= /c, rui e Z>o 



i=l 



i=l 



where the cD^'s denote the fundamental afRne weights of the afRne Lie algebra su(n); see 
e.g. [13] for details. Note that we use the label n instead of for the afRne node. In 
what follows it will be convenient to identify elements in the set with the partitions 
V<n-i,k whose Young diagram fits into a bounding box of height n — 1 and width k. 
Namely, define a bijection — )■ V<n-i,k by setting 

(1.2) A H^. A = (Ai, . . . , An-i) with Ai - A^+i = , 

where mj is the so-called Dynkin label, i.e. the coefficient of the i^^ fundamental weight 
in (11.11) . Vice versa, given a partition A G V<n-i,k we shall denote by A the correspond- 
ing affine weig ht in P+. 

Since the set of dominant integral weights at fixed level k has cardinality \Pjl\ — 
WZNW models are so-called rational conformal field theories, i.e. they have 
a finite set of primary fields from which all other fields can be generated. An important 
ingredient in the description of rational conformal field theories is the concept of fusion: 
in physics terminology one considers the operator product expansion of two primary 
fields. While this can be made mathematically precise in the context of vertex operator 
algebras and the fusion process can be identified with the product in the Grothendieck 
ring of an abelian braided monoidal category in the context of tilting modules of quan- 
tum groups, we will not use this mathematical framework here. 

Consider the free abelian group (with respect to addition) generated by Pj^ and 
introduce the so-called fusion product 

(1.3) ^*A=E-^ir'^' 



by defining the structure constants G Z>o, called fusion coefficients, via the 

celebrated Verlinde formula [191 



(1 4) A/'-^^"' = ^ ^\a^f^^^i>a 



Here denotes the weight corresponding to the empty partition and 5" is the modular 5"- 
matrix describing the modular transformation r — )■ — of affine characters. Among 
other properties it enjoys unitarity Sf^a = S^^ and crossing symmetry S^^ = S^*~, 
where A is the dual weight obtained by taking the complement of A in the (n — 1) x k 
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bounding box and then deleting all n-columns, A* = (Ai, Ai — A^_i, Ai — Xn-2j • ■ ■ , Ai — 
A2). 

For WZNW models the explicit expression for S is known: the Kac-Peterson for- 
mula [H] states S in terms of the Weyl group W and specialises for su(n) to the 
expression 

(1.5) S,, = ^ ^ (_i)^(He-^(-+-(^+'')) , 

v/n(/c + n)^ 

where p is the Weyl vector and A, a denote the finite, non-afRne weights corresponding 
to A, a. From this formula it is by no means obvious that the fusion coefficients ( II ■4p 
are non-negative integers, however they have been identified with certain dimensions or 
multiplicities in various different contexts as e.g. discussed in [9] (for references see loc. 
cit.): dimensions of spaces of conformal blocks of 3-point functions, so-called moduli 
spaces of generalised 6'-functions; outer multiplicities of truncated tensor products of 
tilting modules of quantum groups at roots of unity; Littlewood-Richardson coefficients 
of Hecke algebras at roots of unity; dimensions of local states in restricted-solid-on- 
solid models. In fact, ( II. 3p defines a unital, commutative ring over the integers Z, 
which we shall refer to as the su(n) fusion ring at level k, denoted by J-'n,fc, and to the 
corresponding unital, commutative and associative algebra J-^ J^n,k®z'C as fusion 
or Verlinde algebra. 

This article aims to give a non-technical account of the main findings in [TSj and 
p!6] regarding the iu(n) fusion ring. For proofs the reader is referred to the mentioned 
papers. Sections 2 and 3 are largely a summary of previous results reviewing the defi- 
nition of an integrable statistical mechanics model which generates the fusion ring. It 
is convenient to describe the statistical model and its lattice configurations using non- 
intersecting paths, since this allows for instance a non-technical definition of the transfer 
matrix. However, it needs to be stressed that at the moment the path picture is not 
used to give combinatorial proofs, instead the discussion is algebraic and employs the 
solution to the Yang-Baxter equation given in [16]. However, we present one result. 
Corollary 12.31 which relates the counting of non-intersecting paths on the cylinder to a 
sum over fusion coefficients. We also make contact with the phase model of Bogoliubov, 
Izergin and Kitanine [2] where closely related algebraic structures have been discussed. 
Section 4 states a detailed derivation of the new algorithm to compute fusion coeffi- 
cients. First a review of the Bethe ansatz equations is given by highlighting how they 
are connected to a fusion potential. The latter differs from the familiar fusion potential 
of Gepner [7J and the algorithm therefore yields expressions for fusion coefficients in 
terms of Littlewood-Richardson coefficients which differ from the ones obtained via the 
celebrated Kac- Walton formula [13] [2D]; compare also with the work of Goodman and 
Wenzl [TO]. Section 5 stresses that the Bethe vectors constructed via the quantum in- 
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verse scattering method can be identified with the idempotents of the fusion ring. The 
proof of this result has been given before [15] but their role has not been emphasized. 
The modular S-matrix is the transition matrix from the basis of integrable weights to 
the basis of Bethe vectors and, hence, can be expressed in terms of the generators of 
a Yang-Baxter algebra. We also present a new perspective on the afRne plactic Schur 
polynomials which are defined in terms of the transfer matrix via a determinant for- 
mula: they constitute the set of conserved quantities of the integrable model and hence 
should be seen as the quantum analogue of a spectral curve. In the present model this 
quantum spectral curve coincides with the collection of su{n)k fusion rings for k G Z>o. 
We conclude with some practical applications, recursion formulae for fusion coefficients 
at different level. The last section discusses how the findings summarised here might 
generalise to a wider class of integrable models. 



We start by defining a statistical vertex model which is obtained in the crystal 
limit of UgSu{2); see [16]. Consider a n x (n — 1) square lattice with quasi-periodic 
boundary conditions in the horizontal direction, i.e. a square lattice on a cylinder with 
n — 1 rows. On the edges of the square lattice live statistical variables m G Z>o, which 
we will identify with the Dynkin labels of dominant integrable weights below. To each 
lattice configuration we assign a "Boltzmann weight" in C[;2] [xi, . . . , by taking 

the product over (local) vertex configurations. 

Label the statistical variables a, 6, c, (i G Z>o sitting on the edges of a vertex in the 
i*^ lattice row as shown in Figured! then we assign to it the weight (compare with |16] ) 



It is convenient to describe the vertex configurations in terms of non-intersecting paths, 
so-called oo-friendly walkers; compare with [8l[lTl[T2]. In Figured] 6 walkers are entering 
the vertex from above turning to the right, at which point a contingent of them of size b— 
c chooses to defect. The defectors then join another group of a walkers coming from the 
left. The Boltzmann weight of the i*^ lattice row is then given by 2"" j^K^x^ 
where we have introduced a parameter z which keeps track of how many walkers pass 
the boundary. Figure |3] shows on a simple example that the weight of a single row is 
easily computed by simply counting the number of horizontal edges. 

Given two arbitrary but fixed afhne weights A, G , denote by Tn{p)^ m{P) their 
n-tuples of Dynkin labels in ( II. ip . Denote by the lattice configurations where the 
outer vertical edges at the bottom and top of the cylinder take the values m(/i), m(i>), 
respectively. Figure [2] shows on an example that each lattice configuration corresponds 



§ 2. Fusion coefficients from statistical mechanics 
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Figure 1. Graphical depiction of a vertex configuration with Boltzmann weight ( 12. ip . 
The statistical variables a, 6, c, (i G Z>o obey the constraints a + b — c + d and b > c. 



to k nonintersecting paths some of which are closely bunched together. The corre- 
sponding partition function of the vertex model, i.e. the weighted sum over all lattice 
configurations is defined as 



(2.2) 



E n 

,dij)}erl l<^<n-l 



n 

l<j<n 



'V„, J. . 



with {aij,bij,Cij,dij 



denoting the vertex configuration in the z*'^ row and j^^ col- 
umn. As we will see below the partition function is symmetric in the variables Xi 
and, therefore, can be expanded into a suitable basis in the ring of symmetric func- 
tions C[z][ 

jE5„_i -yy^g choose the basis of Schur functions {sx} where A is a 
partition with length i{X) < n. 

We remind the reader that the Schur function sx{xi, ...,Xn-i) can be defined as 
weighted sum of Young tableaux. Given a partition A, a Young tableau t of shape \t\ = A 
is a filling of the Young diagram with integers in the set { 1 , . . . , ti — 1 } such that the 
numbers are weakly increasing in each row from left to right and are strictly increasing 
in each column from top to bottom. To each tableau we assign the weight vector 
a = (cKi, . . . , ctn-i) where ai is the multiplicity of i occuring in t. The Schur function is 



then given as sx{xi, 



I r 



lit) 



We state an explicit example. 



Example 2.1. Let n = 3 and A = (2, 1). Then the list of possible tableaux t 



reads, 



1 


1 


2 





1 


1 


3 





1 


2 


2 





1 


3 


2 





1 


2 


3 





1 


3 


3 





Thus, the Schur function is the following polynomial 

-5(2,1) = X1X2 + xfxs + X1X2 + 2X1X2X3 + X1X3 -|- X2X3 + X2X3 . 



Expanding the partition function (12. 2p with respect to Schur functions we obtain a 
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relation between the statistical mechanics model defined via ( 12. ip and the fusion algebra 
of the su(n)fc-WZNW model M\. 



Proposition 2.2 (generating function for fusion coefficients). The partition func- 
tion l\2.^) has the expansion 



(2.3) Z^ix,, . . . , a:„_i; z) = ^"-^ff 



A 7 ■ ■ ■ 7 -^n— 1 J 1 



where are the fusion coefficients and the degree d is given 

Given a square s at position (z, j) in the Young diagram of a partition A, recall that 
the hook length is defined as h{s) = + A* — z — j ' + 1. That is, h{s) is the number 
of squares to right in the same row and the number of squares in the column below it 
plus one (for the square itself). The content of the same square is simply defined as 
c(s) = j — i. Denote by r^((i) C the subset of lattice path configurations which 
have a fixed number of 2d outer horizontal edges. The following result on the number 
of possible lattice configurations on the cylinder appears to be new. 

Corollary 2.3 (lattice configurations and fusion coefficients). Specialising to z - 
xi — ■ ■ ■ — Xn-i = 1 in \2.3\) we obtain the identity 

(2.4) |rA(,)|. 5;ArW.^n^^(^. 

A€P+ 



where the sum can he restricted to those weights A for which d — ^■'^^^^^ — H -j_ _ 

Proof. The assertion follows from the well-known formula piSf Chapter I, Section 
3, Example 4], sa(1, • ■ ■ , 1) = nseA('^ — 1 + c{s))/h{s) and the previous expansion ( 12. 3|) 
of the partition function. □ 

Remark. While in [H [TTl [12] the Gessel-Viennot method has been used to obtain 
analogous results for different type of boundary conditions, the proof of Proposition 12.21 
rests on the Yang-Baxter equation, see (13.161) below, and the quantum inverse scattering 
method, thus it is of algebraic nature. 



§ 3. Transfer matrix and Yang-Baxter algebras 

We now introduce the row-to-row transfer matrix of the vertex model (12. ip as the 
partition function of a single lattice row and then identify it below as generating function 
of certain polynomials in a noncommutative alphabet. 
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m(jl) 2 2 10 
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Figure 2. An example of a lattice configuration in terms of non-intersecting paths for 
n = k = 5 and = (5, 3, 1), = (4, 2, 1). 
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1 















Figure 3. An example for n = 5 and /c = 7 of a row configuration and its "statistical 
weight" , which is obtained by counting the horizontal edges of paths. Outer edges are 
identified and for each a power of z is added. 
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Definition 3.1 (transfer matrix). Given any two n-tuples m = (mi,...,m^) 
and m' = {m'l, . . . , m^) G Z>q, tlie transfer matrix Q{xi) of tlie i^^ row is defined via 
ttie elements 

(3.1) Q{x.)m,m':= E ^^-H|^^#ofho..o„taledges^ 

allowed row 
configurations 

where the factor 1/2 in the power of the variable z takes into account that the outer 
horizontal edges need to be identified, since we are on the cylinder. 

As it is common in the discussion of vertex models we wish to identify the transfer 
matrix as an endomorphism of a vector space. For this purpose we now interpret the 
statistical variables at the lattice edges as labels of basis vectors in the vector space 
A4 = ®m^z>o '^^m- Then a row configuration in the lattice, i.e. an assignment of 
statistical variables m = (mi, . . . ,m„) along one row of vertical edges, fixes a vector 

® ■ ■ ■ ® G A^®"". Henceforth, we identify the tensor product A^®"- with CP""", 
the complex linear span of all the integral dominant weights of the affine Lie algebra 
iu(n), P+ := I A = J2i=i ^i'^i "^i ^ ^>o| 1 via the map m i— ^ Xir=i ^i'^i- That is, 
we interpret the statistical variables m = (mi, . . . , m„) in one row of our vertex model 
as Dynkin labels of an affine weight in P+. For convenience we will sometimes denote 
this n-tuple m and the associated vector in Ai®^ by the same symbol. By construction 
the row-to-row transfer matrix and the partition function are then related via 

(3.2) ^/^(a^i, • • • , Xn-i, z) = (m(i>), Q(x^_i) ■ ■ • Q{xi)m{fi)) , 

where we have introduced the inner product (m, m') := Y[i=i ^mi.m'^ which we assume 
to be antilinear in the first factor. Thus, the transfer matrix (13. Ij) can be interpreted as 
discrete time evolution operator which successively generates the paths on the cylindric 
square lattice. Note that for any pair of configurations m,m' G A^®"^ only a finite 
number of the terms making up the matrix element {m,Q{xi)m') is non-zero. The 
operator Q G EndAi®^ is therefore well-defined. We now reformulate the transfer 
matrix in terms of a set of more elementary, local operators which respectively increase 
and decrease a single Dynkin label m^ only. 

§ 3.1. The local affine plactic algebra 

For i = 1, . . . ,n define maps ^i,^* G End(A^®"') by setting 

(3.3) (/?*m = (mi, . . . , mi -Fl, . . . ,m„) 
and 



(3.4) ^^m 



(mi, . . . , mi - 1, . . . , mn),mi > 
0, else 
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In addition let Ni-m = rrii m for all 1 < z < n. These maps can be identified with the 
Chevalley generators of the UqS[{2) Verma module in the crystal limit; see [16]. They 
have first appeared in the context of the phase model; see [2] and references therein. In 
[15] the following statement has been proven by constructing an explicit basis for the 
phase algebra $. 

Proposition 3.2 (phase algebra). The fi^f* and Ni generate a subalgebra $ 
o/End(A^®") which can be realized as the algebra $ with the following generators and 
relations for I < i,j < n: 



(3.5) ip.ip^ = ip^ip,, ip^ip* = ip*ipl N,Nj = NjN, 

(3.6) Nupj - ipjNi = -5ijipi, Nup* - ip*Ni = Sijcp*, 
(3-7) ipi^* = 1, <Pi^* = if]^, if i ^ J, 

(3.8) iV,(l - ip^ip^) = (1 - ^,VJiV, = . 



Note that with respect to the scalar product introduced above we have {flm, m') ~ 
(m, tpiTn') for any m, m' G A^®". 

Definition 3.3 (local affine plactic algebra). Let PI = Pl(^) be the free algebra 
generated by the elements of ^ = {ai, 02, ... a^} modulo the relations 

(3.9) aiOj — ttjai = 0, if |z — j| 7^ 1 mod n, 

(3.10) a^-i-ia^ = Oiai^iOi, af_^iai = Oi^iaiOi^i, 

where (I3.1UI) are the plactic relations on the circle, i.e. all indices are defined modulo 
n. Denote by Plgn = Plfin(-^') the local finite plactic algebra generated from A' = 
{ai,a2, . ..an-i}. 

We recall the following result from p!5] Prop 5.8]: 

Proposition 3.4. There is a homomorphism of algebras Plfin — > $ such that 

(3.11) aj f-^ </j J- , j = l,...,n-l, 

where the representation of the phase algebra $ given by h3. 5]) and \3-4^ lifts to a repre- 
sentation of the local plactic algebra Plfin- Mapping = a„ to ziplip^ it lifts in addition 
to a representation of PI on A4[z]^^ with A4[z] = C[z] ®c o,i^d z an indeterminate. 
Both representations are faithful. 

Remark. The finite plactic algebra first appeared in the context of the Robinson- 
Schensted correspondence: given a word in a noncommutative alphabet it can be 
mapped onto a pair of Young tableaux, usually called (P, Q), by using the bumping 
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algorithm. Q is the recording tableau encoding the sequence of bumping processes; see 
e.g. [6j for an explanation. Lascoux and Schiitzenberger showed that identifying words 
which only differ in their recording tableaux is equivalent to a set of identities of which 
(13.101) are special cases. The local finite plactic algebra was first considered by Fomin 
and Greene in [5]. We recover their case when specialising to 2 = 0. 

Remark. Note that the action of the affine plactic algebra is blockdiagonal with 
respect to the decomposition CP+ = 0^>qCP^. In fact each subspace can be repre- 
sented as a directed coloured graph where the elements in are the vertices and a 
directed edge of colour i between two vertices, fi A, is introduced if A = a^/i. This 
yields the Kirillov-Reshetikhin crystal graph Bi^k of type A. Setting ^ = all edges 
related to the affine generator are removed from the graph and we obtain the su(n) 
crystal graph of highest weight kuji, where ui is the first fundamental weight. 

§ 3.2. Yang-Baxter algebras 

For r, s G Z>o arbitrary but fixed and u a formal, invertible variable define Q{u) : 
M{u) ^Miu)^"^ M{u) ^Miu)®"^ by setting Q{u)vr ®m := Zls>o^s (8) Q^^^(w)m 
where 

(3.12) Q,^M :=E^'''^'^(^i)''«i'---<-"iV^ 

e 

with the sum running over all compositions e = (ei, . . . , En-i)- Despite the sums in the 
definition of Q{u) being infinite, only a finite number of terms survive when acting on 
a vector in Ai{u)^^~^^ , thus the operator is well-defined. In fact, we have the following 



Lemma 3.5. Let Q be the transfer matrix defined in ( 15'. Then 

(3.13) Q(x.) = J2 ^''QrA^^) = $^(^x,)'-(^D'^Qo,o(^.)^;, 

r>0 r>0 

in other words for z = 1 the transfer matrix is the formal trace of the matrix liS.l^) . 
Define another operator lZ{u/v) : M.{u) ® M.{v) M.[u) ® J^iy) via the relation 

(3.14) n{u) Va®Vb= 



,d>0 



setting 



w", c = b, d = a 

(3.15) ^c.dC") = { -u),d = a + b- c, b>c 

0, else 
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Figure 4. The (Boltzmann) weights of local vertex configuration for the auxiliary matrix 
T = A + zD. Note that in contrast to the Boltzmann weights (12. II) the walkers can now 
propagate horizontally but only one is allowed on a horizontal edge at a time. 



Proposition 3.6. The TZ-matrix I!i3.14\ ) o.'^^d the monodromy matrix ^3.12^ obey 
the Yang-Baxter equation, 



(3.16) 



1li2{u/v)Qi{u)Q2{v) = Q2{v)Qi{u)TZi2{u/v) . 



The definition of the monodromy matrix ( I3.12p can be algebraically motivated 
by taking a special limit of the intertwiner of a Uq{5l{2)) Verma module; see the dis- 
cussion in [16]. Replacing this Verma module with the (two-dimensional) fundamen- 
tal representation one obtains in the analogous limit now a 2 x 2 monodromy ma- 
trix, T G End(C^ ® A^(m)®"), where the ordering of the noncommutative alphabet is 
reversed(r, s = 0, 1), 



(3.17) 



" 'f'n"n-l 



A{u) B{u) 
C{u) D{u) 



Note that the sum now only runs over compositions e whose parts are or 1. This 
second monodromy matrix coincides with the matrix introduced by Bogoliubov, Izergin 
and Kitanine in the context of the so-called phase model [2]. 



u 



Proposition 3.7 ([2J). The 2x2 matrix ( [g.j7| j with entries in End(A1( 
obeys the Yang-Baxter equation Ri2{u/ v)Ti{u)T2{v) = T2{v)Ti{u)Ri2(u/v) with 

(■ 

(3.18) R{u) = 



u 
u-1 

















u 
u-1 








1 

u-1 


1 














u 
u-1 
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Figure 5. An example of a row configuration for the auxiliary matrix T for n = 6 and 
k = 7. The statistical weight is again obtained by counting the horizontal edges of 
paths, where outer edges are once more identified each contributing a power of z. 



The matrix elements of (I3.17P obey certain commutation relations with the matrix 
elements (I3.12p , which again can be encoded in yet another solution to the Yang-Baxter 
equation [16]. 



Proposition 3.8. The monodromy matrices li3.1S\) and p.jTp satisfy the equa- 
tion, 

(3.19) L[^{u/v)Ti{u)Q^{v) = Q^{v)Ti{u)L\^{u/v), 

where the endomorphism L' G End(C^ ® M.{u)'^^) is defined in terms the matrix ele- 
ments 

{1 + M, a, 6, c, (i = 

w", d = a + 6-c, 6 >c, a, c = 0,1 . 
0, else 

The hallmark of an exactly solvable or integrable model in statistical mechanics 
is that its transfer matrix commutes with itself for arbitrary values of the spectral 
parameter which here is identified with the variables Xi in each lattice row. In a physical 
application the row variables {xi, . . . , x^} would be evaluated in the interval [0, 1] ^"^ 
such that the Boltzmann weights ( 12. ip can be interpreted as proper probabilities. The 
transfer matrices for any other, possibly complex values, of the XiS would be seen as a 
"symmetry" of the system. Generalising the notion of Liouville integrability in classical 
mechanics, such a statistical model is called integrable. One important consequence of 
the Yang-Baxter equations stated above is that they imply integrability of the vertex 
model (1^ . 

Corollary 3.9 (Integrability). SetT = A-\-zD then we have, among others, the 
commutation relations 



(3.21) [Q{u),Q{v)] = [nu),T{v)] = [T(w),g(t;)] =0, 
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where Q is the transfer matrix liS. Moreover, the following relation holds true 

Qr,MB{v) = (l + ^) 6r,oB{v)Qo,o{u) 
(3.22) +Q^_,^^iu)Div) - Q,^,+,iu)A{v) + Q,_^^,+^{u)C{v) 

Note in particular that the first relation in ( I3.2ip entails via (13. 21) that the partition 
function (12 .21) is symmetric in the variables claimed earlier. 

§ 3.3. AfRne plactic elementary and complete symmetric polynomials 

To keep this article self-contained and motivate the definition of the affine plactic 
polynomials below, we review some basic facts about symmetric functions; see [18j 
for details. Let {yi, . . . ,ye} be a set of commuting variables for some finite i > 0. 
Recall that the ring of symmetric functions C[?/i, . . .,yi\^^ is generated by either the 
elementary or complete symmetric functions denoted by the letters and hr with r G 
Z>o, respectively. Both sets of functions can be introduced via the following generating 
functions, 



(3.23) JJ(l + 2/iw)=^e^(yi,...,y£)w'' 

1=1 r>0 

I 

(3.24) - ^^^)"' = 5Z ^-(^1' • • • ' 2/^)"^ 

i=\ r>0 

respectively. Note that the first sum is finite, i.e. er{yi, ■ ■ , ye) = for r > i, while the 
second one is infinite. Explicitly, the elementary and complete symmetric functions are 
given by the expressions 

(3.25) er{yi,...,ye)= ^ Vh-'-yir- 

l<ii<---<ir<i 

(3.26) hr{yi,...,ye)= ^ yir'-'Vii- 

l<ii<---<ir<e 

Prom the generating functions it is immediate to deduce that the e's and /I's satisfy 
recursion relations. Namely, one has the identities 

(3.27) er{yi,...,yi) = er{yi, . . . ,yi-i) + yiCr-iiyi, . . . ,ye-i), 

(3.28) hr{yi,...,yt)=hr{yi,..., ye-i) + yehr-iiyi, ■ ■ ■ ^yt) ■ 

As mentioned above one has C[yi, . . . , yif''^ = C[ei, . . . , e^] = C[/ii, . . . ,hi\ and, thus, 
both sets of functions must be polynomials of each other. In fact, multiplying both 
generating functions (with u replaced by —u in the first one) yields the equations 
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J2r=a+b(^^)'^^a^b = for all T > which can be solved either for the e's or the 
/I's to yield the well-known Jacobi-Trudi formulae 

(3.29) Kivi, ...,yi) = det(ei_i+j(yi, . . . , ye))i<i,j<r: 

(3.30) erivi, ...,yi) = det(/ii_i+j(yi, . . . ,yi))i<i,j<r ■ 

We will now generalise these functions by replacing the commuting variables {yi, . . . , 
yi} with the noncommutative variables {ai, . . . , a^}, i.e. the generators of the affine 
plactic algebra. We wish to identify the auxiliary matrix T defined in Corollary 13. 9[ 
the trace of (13.171) . and the transfer matrix Q defined either through ( 13. ip or ( I3.13P as 
noncommutative analogues of the generating functions ( 13.230 and ( I3.24p . respectively. 
When trying to generalise ( I3.25P and (I3.26P to a noncommutative alphabet one needs 
to specify an ordering of the variables. The auxiliary and transfer matrix prescribe 
such an ordering of the noncommutative variables {ai} which turns out to be consistent 
and has the additional desirable property that the resulting affine plactic elementary 
and complete symmetric polynomials form a commutative subalgebra due to the inte- 
grability condition (13.211) . Because this ordering is cyclic, it is easier to split the affine 
generator a^ into its constituents z, Lp\ and ip^, compare with Proposition 13. 4[ and write 
the variables in descending or ascending order as we have already done in (I3.17P and 
(l3J2il . 

Proposition 3.10 (generating functions). The transfer matrix \3. 1\] can he in- 
terpreted as the generating function 

(3.31) Q{xi) = z'QrA^i) = Yl ^"hriA) 

r>Q r>0 

of the affine plactic complete symmetric polynomials 

(3.32) K{A) :=E^'''(^i)'"«i'---<-iVn° • 

ehr 

The auxiliary matrix defined in Corollary \3.9\. 

(3.33) T(xO = A(x,) + zD{x,) = J] x[e,(^), 

r>0 

on the other hand, yields the generating function of the affine plactic elementary sym- 
metric polynomials, 

(3.34) eM)--= E 

ehr 
e«=0,l 
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§ 3.4. Combinatorial description of the Yang-Baxter algebra 

Note that when setting the quasi-periodicity parameter z to zero, i.e. enforcing open 
boundary conditions, one obtains the finite plactic polynomials of [5]. For instance, we 
have that the matrix element A in ( I3.17P is the generating function of the finite plactic 
elementary polynomials, 

(3.35) ;z = 0: = (1 + a^-iw) ■ ■ ■ (1 + aiw) = ^ w'^e^(^') ■ 

r>0 

Lemma 3.11. The action of the operator er{A') = er(ai, . . . , a„_i) on a parti- 
tion li G V<n-i,k produces a sum over all partitions X in the n x k bounding box such 
that the skew diagram X/n is a vertical strip of length r and Xi = Hi, 

(3.36) eriA')^^= J2 '^A^M.A'. 

Here X' G V<n-i,k is the partition obtained from X by deleting all columns of height n. 

Proof. Using the graphical depiction of the Boltzamnn weights in Figure H] it 
follows that for z = only row configurations are allowed which do not have an occupied 
outer horizontal edge. Hence, mi{fi) > mi{i/) which entails that we must have fii = 
ui according to (11.21) . Moreover, we can deduce from the Boltzmann weights that 
mi{u) — mi{fi) = or 1 for i > 1, hence A//i must be a horizontal r-strip. □ 

In order to describe the action of the remaining matrix elements in (13.171) note that 
in terms of partitions the map <f>l adds a column of height one and increases the width 
of the bounding box, the level k, by one. The map (p^ simply decreases the width of 
the bounding box if Ai < /c, otherwise it sends A to zero. Thus, the following formulae 

(3.37) B{u) = u A{u)^l, C{u) = ^^A{u), D{u) = u^^A{u)^l , 

which can be easily checked from (13.171) . allow one to perform computations with the 
Yang-Baxter algebra purely in terms of Young diagrams and their bounding boxes. 

Example 3.12. Let n = h and choose n = (2,2,1,0) with m = (0,1,1,0,1), 
that is k = 3. Then we have for r = 3 only a single term. 
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In contrast the action of the ajfine plactic polynomial yields the sum 



+ z 



+ z 









+ z 

































+ z 



By converting each partition in the above sum to the corresponding compositions one 
verifies that each of the last four terms on the right hand side is generated by a monomial 
in A which contains the affine generator an- For instance, for the second term one finds 



z(pr^a2ai(pi 



m=(0,l,2,0,0) 



m=(0,l,l,0,l) 



where we have used (13 .71) to rewrite the respective term in es{A) as word in the affine 
plactic generators. This is always possible for r < n. 



In light of (13.351) and the last example the definition of the auxiliary matrix ( I3.33P 
can be seen as the noncommutative analogue of ( I3.27p . since after expanding with 
respect to the variable u one arrives at the identity 



(3.38) 



er{A) = er{A') + z(p^er-i{A')(pl . 



Cylic ordering. As already alluded to in the last example comparison with the com- 
mutative case ( I3.27P is made easier by realising that for r < n the terms in the second 
summand can always be rearranged in cyclic order. To expose the general structure 
more clearly consider another example provided by the row configuration depicted in 
Figure [5] for n = 6 and k = 7. The latter corresponds via ( I3.34p to the monomial 

zipQa5a4a2ai(pi = za2aiLpQa^a/^Lp\ = 0201050504, 

where we have used once more ( 13. 7p . The general case is now clear: monomials in the 
tti which do neither contain oi or o^, or only one of these generators, are always written 
in descending order from left to right. If both, oi and o„, occur in the same monomial 
write the maximal string of form o^o^-i ■ ■ ■ 02O1 to the left of the remaining letters which 
should also be in descending order starting with o^ (although the indices might now 
"jump" by more than one). It follows from the definition ( I3.34p that for r = n we have 
en{A) = z-1. 

We now specialise to the finite plactic complete symmetric polynomials, 
(3.39) z = 0: Q{u) = Qo,o(«) = Yl ^^^M^') ■ 

r>0 
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From the definition ( I3.32p one easily computes the noncommutative analogue of the 
recursion relation ( I3.28p . 

(3.40) hr{A) = KiA') + ziplhr-i{A)ip^ . 
Lemma 3.13. Let r < k and ^ G V<n-i,k then 

(3.41) K{A')ii= '^Ai-^i^'- 

For r > k we have hr{A')\^p+ = 0. 

Proof. The assertion follows from a similar line of argument as before, but this 
time one uses the Boltzmann weights depicted in Figure [H Since z = row config- 
urations with outer edges are prohibited, whence mi(//) > mi(z^). In contrast to the 
previous case (I3.36P a friendly walker now cannot propagate horizontally, however sev- 
eral are allowed at the same time on the horizontal edges. Thus, we obtain a horizontal 
instead of a vertical strip. The last identity is also clear from the graphical depiction of 
allowed row configurations: the number of occupied horizontal edges cannot exceed the 
number of incoming walkers. □ 

Note that the affine plactic complete symmetric polynomials can only be rewritten 
in (reverse) cyclic order for r < n using the same commutation relations of the phase 
algebra as before. For r > n the cyclic ordering ceases to be well-defined and one has 
to resort to f l3.32p . 

Finally, we generalise the last identities from the commutative case, the Jacobi- 
Trudi formulae (13.290 and fl3.30p . which are subject of the next proposition [16]. 

Proposition 3.14 (operator functional equation). The generating functions liS. 31\) 
and lis. 33\) satisfy the operator functional relation, 

(3.42) T{-u)Qiu) = 1 + z{-ir Yl «'+"^fc(^)7rfc , 

fc>0 

where tt^ : CP~^ — )■ CP^^ is the projector onto the subspace spanned by the weights at level 
k. In particular, for r < n + k the familiar determinant relations from the commutative 
case also hold for the noncommutative elementary and complete symmetric polynomials, 

(3.43) hr{A) = det{ei-i+j{A))i<ij<r, er{A) = det{hi_i+j{A))i<ij<r , 

where the determinants are well defined due to /i3.21\) . 

Setting once more z = we recover the relation expected from the commutative 
case. The additional terms have their origin in the quasi-periodic boundary conditions 
and we explain their origin on an example which will elucidate the general formula. 



18 



Christian Korff 



















Figure 6. The two possible vacuum row configurations of the auxihary matrix T. Note 
that the transfer matrix Q only allows for the left one. 



Q 

D 



10 



10 



Figure 7. Depicted is the row configuration resulting from the successive action of the 
transfer matrix Q and the operator D on a state of level k = 1 and which is responsible 
for the additional term on the right hand side of the functional equation (I3.42p . 
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Example 3.15. Set n > 2 and consider first the case of level k = 0. There is only 
one state, the "pseudo-vacuum" O = (0,0,.. .,0), and trivially we have Q{u)Q = O. 
Acting with T{—u) on we obtain two contributions shown in Figure El Thus, the 
additional term in f l3.42p with Hq^A) = 1 is due to a "vacuum mode", a path which 
winds around the cylinder. 

Let us now set /c = 1 and consider the state Ui = (0, . . . , 0, 1, 0, . . . 0) with i < n. 

i 

Then Q{u)Cji = Ui + wwi+i as can be deduced graphically from ( 13. ip . In order to 
understand the origin of the additional term in (I3.42p which includes the factor z, it 
suffices to look at the contribution from the D operator, since T = A + zD and A, D do 
not contain z. The action of D can be easily computed using Figure H] and remembering 
that only row configurations contribute where the outer edges are occupied. One finds 
that all row configurations cancel except for one term which is depicted in Figure [71 
again the additional term in ( I3.42p corresponds to a diagram where one path winds 
around the cylinder. 

§ 3.5. Quantum Hamiltonian and particle picture 

So far we discussed a statistical mechanics model whose transfer matrix can be iden- 
tified with the generating function of the affine plactic complete symmetric polynomials. 
Similarly, we can link the generating function of the affine plactic elementary symmetric 
polynomials, the auxiliary matrix T, to a physical model in quantum mechancis. Inter- 
pret the Dynkin labels rrii as occupation numbers of a site z of a circular lattice - the 
Dynkin diagram of su(n) in Figure [8]- and the maps ( 13. 3 p and ( 13. 4p as particle creation 
and annihilation operators, respectively. Introducing the quantum Hamiltonian 

(3.44) H = = -- + z'^j'^j+i) 

3 = 1 

with = z~'^(fii, ^n+i = z^i^ '^^^ the conserved charges Hf" = —{er{A) ± 

en—r{,A)^ /I which because of (|3.2ip are in involution, [H,Hf] = [H^,H^] = 0, yields 
an alternative physical interpretation of the combinatorial structures described here. 
This quantum system is known as phase model, see [2J and references therein. 

§ 4. Bethe ansatz equations and the fusion potential 

Within the framework of exactly solvable models the next step is to construct the 
eigenvectors of the transfer matrix Q. Instead it is simpler to consider the eigenvalue 
problem of the auxiliary matrix T, since it follows from the functional relation (13.420 and 
the determinant formulae (13.430 that the eigenvectors of Q coincide with the eigenvectors 
of T. The advantage of this approach is that the eigenvectors of T can be computed via 
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Figure 8. The fusion ring can be interpreted as a discrete model in quantum mechanics. 
Each Dynkin label rui specifies a number of particles sitting at the i*^ node of the 
Dynkin diagram. Application of the generator of the affine plactic algebra, shifts one 
particle from site i to site i + 1, as indicated in the picture for z = 1. 

the algebraic Bethe ansatz or quantum inverse scattering method (see e.g. j3] for a text 
book and references therein) using the commutation relations of the A, B, C, D algebra 
in (I3.17p . which are drastically simpler than the commutation relations of the algebra 
generated by the matrix elements ( I3.12p . Starting point is a particular assumption on 
the algebraic form of the eigenvectors: define an (off-shell) Bethe vector at level k > 
to be 

(4.1) 6(xi, . . . , Xfc) := S(xr') ■ ■ ■ B{x-^)Cl, 

where O is again the unique vector corresponding to the composition (0, 0, . . . , 0). The 
requirement that (14. ip is an eigenvector of T leads via the commutation relations of the 
Yang-Baxter algebra contained in ( I3.18P and a standard computation - which we omit 
- to the Bethe ansatz equations |2j [15\ 

k 

(4.2) = • ■ ■ = x^+'' = II X, . 

i=l 

We now discuss how the Bethe ansatz equations lead to a combinatorial computa- 
tion of fusion coefficients. We wish to emphasize that this is possible without solving 
(14. 2 p first. For this reason we postpone the discussion of their solutions to the next 
section, however, we already mention that in order to solve them one needs to assume 
that ^ii/"- exist. 

The Bethe ansatz equations are k polynomial equations and, thus, describe an affine 
variety ^ C C[2;^"]'^. Recall that given a field K an affine variety is usually defined 
in terms of a set of polynomials /i, . . . , /jt G . . . , Xk] by setting V(/i, . . . , fk) '■ = 

|f G K'^ : /i(f 1, . . . , Vk) = 0, z = 1, . . . , /c}. Note that due to the commutation relations 
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B[x)B{y) = B{y)B{x), which foUow from the Yang-Baxter equation (13. 181) . we can iden- 
tify solutions of ( 14. 2p under permutations (xi, . . . , x^) ~ {xyj^, . . . , x^,^) for aU w G Sfc, 
since they give rise to the same eigenvector (14. 11) . Denote by Vn,fc the variety obtained 
under this identification, then we can think of Nn,k as being defined by k elements in 
the ring of symmetric polynomials C[z^"]^[xi, . . . ,a;fc]^*' = Cfz^^j'^fei, . . . ,efc], where 
the Ci's are the elementary symmetric functions. 

Lemma 4.1. Let hr = det(ei_i_|_j)i<ij<r, the complete symmetric functions, 
then the ajfine variety Yn,k defined by the Bethe ansatz equations (up to permutations 
of the solutions) is given 6j0 

(4.3) Nn,k = V(/ln - Z, hn+l, . . . , hn+k-l) ■ 

The proof of this statement can be found in jl5| Lemma 6.3] and simply uses 
the recursive relation ( I3.28P for complete symmetric polynomials. Next we assign to 
the solutions of the Bethe ansatz equations ( 14. 2 p an ideal in the ring of symmetric 
functions. Recall the definition of the nuUstellen or vanishing ideal of an affine variety 
V, X(V) := {/ e C(z)[ei,...,efc] : f{v) = 0, v e V}. Then we have the following 
statement [151 Proof of Theorem 6.20]: 

Proposition 4.2. LetX(Vn,k) be the nuUstellen ideal of the Bethe ansatz variety 
^73\ ), then 

(4.4) {hn - Z, hn+l, . . . , hn+k-l) = ^'^n,k) ■ 

This proposition is proved via Hilbert's NuUstellensatz which asserts that given an 
ideal / in a polynomial ring one has X(V(/)) = yj, where 77 = {/ : /"^ G / for some 
m G Z>o} is the radical of /. In the present case where / = {hn — 1, /i^+i, • ■ • , hn+k-i) 
one shows that / = y/l and thus the assertion follows from the previous lemma. 

Remark. For z = 1 the ideal fl4.4p can be encoded into a fusion potential pk+n = 
X]i=i 3^?"'"'^ 5 the {k + power sum, noting that 

1 dpk+n , 



(-l)'-^/ifc+„_, = 0, r = l,2,...,A;. 



/c + n dcr 

This is very similar to the fusion potential introduced by Gepner. The difference lies in 
the constraints imposed on the variables: Gepner's fusion potential [3, Equation (2.31)] 



(4.6) y„,, = _L.5^,.".^ 



^In [15j an additional relation has been stated to facilitate the comparison with the small quantum 
cohomology ring of the Grassmannian, = zi^—\y^^^ey.. This last relation follows from /i^—i; = 

^n + l = ■ • • = Yin = by exploiting the definition the hook Schur polynomial S(„|fe_i) = 

e^hn = ze^; see (|4.10p below. 
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is defined in terms of n variables {yi, . . . , y^} subject to tlie constraints e„(y) = yi ■ ■ ■ Vn = 
1 and 

(4.6) = {-iy+'hr,+k-^ = 0, z = 1, . . . , n - 1 . 

Tliese constraints can be sliown to be equivalent witli tfie following set of equations, 

(4.7) = ■ ■ ■ = Vn^' = hk+n{yi, ...,yn) = (- 1)"" ^ /^fc (l/l , . . . , ^n), 

which look very similar to the Bethe ansatz equations (14.21) . However, the corresponding 
construction of the Bethe vector in terms of a Yang-Baxter algebra is currently missing. 

The importance of the ideal ( 14. 4p derived form the Bethe ansatz equations lies in 
the fact that it yields a presentation of the Verlinde or fusion algebra in the ring of 
symmetric functions p!5l Theorem 6.20]. 

Theorem 4.3 (Korff-Stroppel). Set z = 1. Then the map 9 A t— ?■ + 

I(^n,k)j where I{Yn,k) is the ideal in ^.4\) provides an algebra isomorphism, 

(4.8) ®z C = C[ei, . . . , efc]/X(V„,fc) . 

In contrast Gepner's fusion potential is equivalent to a different presentation; c.f. 
[TOl p247, result (4)] and [3 Equation (2.36)]. 

Theorem 4.4 (Gepner, Goodman-Wenzl). The map 9 A h-> SA + /n,A;, where 
k \Cn — 1, hn+iT ■ ■ T hn+k-i) is the ideal resulting from ( [^.g| ) also provides an 
isomorphism, 

(4.9) J-„,fc®zC = C[ei,...,en]//n,fc . 

Proof. For the sake of completeness we briefly outline a proof of ( 14. 9p by showing 
that the ideal (14.61) following from Gepner's fusion potential is identical with the ideal 
used by Goodman and Wenzl in [10] (with the extra condition Cn = 1) who proved 
that the fusion ring is isomorphic to a certain representation of the Hecke algebra at a 
primitive (n + root of unity. 

Let Jn^k be the ideal generated from e„ — 1 and the Schur polynomials of the 
form S(Ai,A2,...,A„) with Ai — A„ = 1; c.f. [TOl p247, result (4)]. Both ideals can be 
shown to be radical along similar lines as it is discussed for (14.41) in Proof of 
Theorem 6.20, Claim 1]. Hence, employing the NuUstellensatz twice, X(V(/„ fc)) = ^ 
and X(V{Jn,k)) = Jn,kj it suffices to prove the two inclusions Jn^k ^ ^{^{In,k)) and 
In,k ^ I{^{Jn,k))- FoT this purpose we recall the definition of hook Schur polynomials 
[m Chapter I, Section 3, Example 9] 

(4.10) S^a\b) ■= S^a+1,1'') = ha+lCb - ha+2eb-l H h {-ifka+b+l ■ 
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Using the Frobenius notation A = (cti, . . . , ar|/?i, • • • , /^r) ^ partition A, where ctj 
and /?j are the lengths of the horizontal and vertical part of a hook centered at the i^^ 
box in the diagonal of A, one has 



(4.11) 



We first show that J C X(V{In,k)) = Exploiting e„ = 1 and the Fieri rule 

enS(Ai,A2,...,A„) = S(Ai+i,A2+i,...,A„+i) we can restrict ourselves to Schur polynomials 
of the form 5(^+1 >2,...,A„_i,o)- From the definition of In,k it then follows that 



\k\b) 



hk+iGb — hk+2eb+i + ■ ■ ■ + {—l)^hk+b+i — 



for all b = 0,1,2, ... ,n — 2 and, thus, we can conclude with the help of (I4.1ip that 
■S(fc+i,A2,...,A„_i,o) = for aU /c + 1 > A2 > ■ • • > A^_i as required. 

The converse inclusion ^ C X(V(J^^fc)) is easily derived along similar lines: 

= S(A;|o) = S(A;+1,0,...,0) = hk+1, 

= ^ik\l) = ^fc+lCl — hk+2 = ~hk+2, 



= S(A;|n-2) = ^fc+ien-2 — ■ ■ • + ( — 1)" '^hk+n-1 = ( — 1)" • 

This shows that the ideal of Goodman and Wenzel coincides with the ideal of Gepner 
and using [101 p247, result (4)] we arrive at f lOll . □ 



Remark. The two isomorphisms (14. 8p and (14.91) lead to different expressions for 
the fusion coefficients in terms of Littlewood-Richardson coefficients. We will discuss 
the case (14.81) below. Goodman and Wenzl used their presentation (14.91) to derive the 
Kac- Walton formula [13] [20] (compare with [10, p247, result (6)]), 



(4.12) 



wew 



where W denotes the affine Weyl group, e{w) the signature oi w, w ■ i> = w{i> + p) — p is 
the shifted Weyl group action with p = J2i '^i being the affine Weyl vector and {w ■ P)' 
is the partition obtained under the bijection (11.21) . 

Example 4.5. Set ti = 3 and /c = 4 and consider the affine weights A = cDq + 
2a)i + Cj2, pi = Cjq + Cji + 2Cj2 in . The corresponding partitions under (11.21) are 
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A = (3,1) and fx = (3,2). Employing the Littlewood-Richardson rule (see e.g. [6]) 
yields the partitions 



(4.13) p = (6, 3, 0), (6, 2, 1), (5, 4, 0), (5, 3, 1), (5, 3, 1), (5, 2, 2), (4, 4, 1), 

(4, 3, 2), (4, 3, 2), (3, 3, 3), (5, 2, 1, 1), (4, 3, 1, 1), (4, 2, 2, 1), (3, 3, 2, 1) . 

Discarding all partitions of length > n and removing all n-columns from the partitions 
u we obtain the following su(n) tensor product decomposition 

A ® ^ = (6, 3) © (5, 1) © (5, 4) © 2(4, 2) © (3, 0) © (3, 3) © 2(2, 1) © (0, 0) . 

Here we have identified highest weight modules with the corresponding partitions. We 
wish to consider the fusion coefficient of the affine weight i> = 2cDi + 2a;2 with partition 
u = (4, 2). From iKWl we find 



(4.14) 



Xjj, 



(4,2) _ (6,3) 



1 = 1, 



since sq ■ (6, 3) = sq ■ (— 2a;o + 3cDi + 30)2) = with sq denoting the affine Weyl reflection. 
In fact, the entire fusion product expansion is computed to 



(4.15) 



+ 



+ 









+ 2 

















+ 



§4.1. Algorithm to compute fusion coefficients 

Based on the presentation (14. Sp derived from the Bethe ansatz equations we now for- 
mulate an alternative algorithm how to compute fusion coefficients in terms of Littlewood- 
Richardson numbers. 

1. Compute the expansion Sx^Sf^t = J2pt c^^Spt via the Littlewood-Richardson rule ; 

note that c^^ = c^t^t Discard all terms for which the partition has length > 
k. 

2. For each of the remaining terms with p\ > n make the replacement Spt = S(^pt^^__^pt^^pt^_ny 
Whenever (p|, . . . , p\, p\ — n) is not a partition use the straightening rules p!8] 

S{...,a,b,...) = -S(...,b-l,a+l,...) and S(. = 

for Schur polynomials to rewrite S(^pt^^^^^^pt^^pt^_^') as Sj^t with i^* a partition. 

3. Collecting terms for each ly one obtains the fusion coefficient J\f['l^'^ . 
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Example 4.6. As in example ( 14 .Sp set n = 3, /c = 4 and consider the parti- 
tions A = (3, 1, 0) and fx = (3, 2, 0). After taking the transpose partitions in (14.131) we 
discard = (2,2,2,1,1,1), (3,2,1,1,1,1), (2,2,2,2,1), (3,2,2,1,1), (3,3,1,1,1) and 
(4, 2, 1, 1, 1) as they have length > k = 4. We are left with three partitions for which 
p\ > 71, namely (4,2,2,1), (4,3,1,1), (4,3,2,0). Employing the above algorithm we 
calculate 

5(4,2,2,1) = ■5(2,2,1,1), 5(4,3,1,1) = S(3, 1,1,1), ^(4,3, 2,0) = ^(3, 2,0,1) = . 

Removing all rows of length n = 3 and collecting terms, one arrives after taking the 
transpose again at the expansion (14.151) . however, the expression for the fusion coef- 
ficients in terms of Littlewood-Richardson numbers are different. For instance, the 
coefficient of = (4, 2, 0) is according to the Kac- Walton formula the difference of two 
Littlewood-Richardson coefficients, see (I4.14p . In contrast, here we find that 



./-(fc)(4,2,0) ^ „(4,3,1,1) ^ . 

•'^(3,1,0)(3,2,0) ^(3, 1,0,0), (3, 2,0,0) 



Thus, the Bethe ansatz equations ( 14.21) provide an alternative algorithm to compute 
fusion coefficients. 



Proof of the algorithm. Consider the ring of symmetric functions C[xi , X2, . . . , 
then s\ = unless i{X) < k. This justifies Step 1 of the algorithm. To deduce Step 
2, assume we are given a partition A with £(A) < k and Ai > n. Then we rewrite the 
Schur function as ITSl 



k 



SA(,Xi,X2,...,Xfc 



weSk l<i<j<k ^' * 



Observing that the equations (14.20 for xi are equivalent to 

x-ll^^=zll^^ 

.X. .X. nr' ^ ■ A A . nr* 

J>1 ^ ^ J>1 ^ ^ 



one derives the identity 



^1 



■xl''Q{x) = zwo- {x^^x^' ■ ■ ■xl'~''e{x)) 



where wq = ak-i ■ ■ ■ cr2cri and is the transposition which permutes Xi and Xj+i. 
Insertion of this identity into the above expression for the Schur function proves that 
sx = ■S(A2,A3,...,Afe,Ai-n) +X(V„,fc). Step 3 then follows from (gS)). □ 

§ 5. Bethe vectors as idempotents 

We now solve the Bethe ansatz equations (14.21) explicitly, which is possible due to 
their simple form, and describe the variety k which consists of a discrete set of points 
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in C[z^/"]'^. For each partition a in V<n-i,k define the following tuple of elements in 

(5.1) P<n-i,fc 3a^x^ = z^C^ (C'^(^'), . . . , C''=^"')), 
where ( = exp(^^) and the exponents are the half-integers, 

(5.2) /(a*) = (^+4-fc,...,^ + al-l) . 

A straightforward computation shows that x'^ solves ( 14. 2 p for any a p!5l Theorem 6.4]. 

Theorem 5.1 (completeness of the Bethe ansatz). Fix k > 0. Then the set of 
vectors := b{xi, . . . , x^)}cr6'P<„„i ^ defined in terms ( [^.ip and Ii5.1\) forms an eigen- 
basis of the transfer l\3. 1\) and auxiliary matrix l{3.33\) in the subspace CP^. In partic- 
ular, let S denote the modular 5u{n)k S -matrix U.5\) . then the Bethe vector b^ has the 
expansion (for simplicity we now label S -matrix elements with partitions) 

(5.3) b.=z-' y: ^'^-^^a. 

Moreover, one has the following eigenvalue equations for the affine plactic polynomials 

(5.4) hr{A)b^ = z^-'^+^ %^ b^ and er{A)b^ = z^^^^^ %^ b^ , 

where (r) and {V) denote the partitions whose Young diagrams consist respectively of a 
single row and a single column of length r. 

Inherent in the last result is the statement that the modular S-matrix can be 
computed in terms of scalar products of the on-shell Bethe vectors b^r and, hence, 
ultimately in terms of the Yang-Baxter algebra generator B via (14.11) . Namely, from 
(15.31) we obtain for 2; = 1 that (60-5 A) = Sx^j/S^a- and {bu.ba-) = |S'0o-|~^. Using that 
500- > we find 

(5.5) z = l: Sx. = ^^^^, b^ = B{xl)---B{xl)n. 



{ba,ba) 



2 



Note in particular, that for a = we obtain the groundstate 60 = X^AeP^ ^^x/ ^ 
of the quantum Hamiltonian (]3.24p . or equivalently the Perron- Frobenius eigenvector 
of the transfer matrix (13. ID . whose components are given by the so-called quantum 
dimensions 



(a,P + A) (a,p + A) 



(5.6) J] C ' -C - ^^^-^ 
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Here the product runs over all positive roots of su(n) and p = ^ Sa>o Weyl 
vector. The expression in terms of Schur functions generalises to the excited states, 
we have in general that Sxcj/S^jj^ = Sxt{x'^) which can be interpreted as characters 
evaluated at special points; see [T5j for details. 

Since we have skipped the algebraic Bethe ansatz computation for T let us verify 
for the simple case k = 1 that the Bethe vector (14.11) is indeed an eigenvector of the 
transfer matrix subject to the Bethe ansatz equations (14.21) and compute the modular 
S-matrix. 

Example 5.2. Take n > 2 and set k = z = 1. Then it follows from ( I3.22p that 
Q{u)B{v)n = (l + ^) B{v)Qo^o{u)n 

+ J2 Qr,r+l{u)mv) - A{v))n + Qr,r+2{U)C{V)(1 
r>0 r>Q 

Exploiting that C{y)Ct = and Qq Q{u)Ct = O, we need to choose the variable v such that 
the second summand on the right hand side vanishes. Observing that {D{v) —A{v))fl = 
[v"^ — 1)0, we arrive at the Bethe ansatz equations (14.21) with xi = v~^. The solutions 
are easily found to be xi{s) = C^^^C^ = with s = 0, 1, . . . , n — 1 and the Bethe 

vector thus reads 

n-l 

h :=S(a:i(s)-i)0 = 5^e'^'^^ {V) . 

r=0 

The modular S-matrix is then easily computed to be 

6 " 



n 



Given the eigenvalue equations ( 15. 4p it is natural to define affine plactic Schur 
polynomials via the familiar Jacobi-Trudi formula (we exploit once more the integrability 
condition ( I3.2ip which guarantees that the determinant is well-defined), 

(5.7) sxiA) := det{hxi-i+jiA))i<ij<n-i ■ 

It is then not difficult to show that the affine plactic Schur polynomials satisfy the eigen- 
value equation sx{A)bcr = {Sxa/ S0a)b(T which leads to the next result [151 Proposition 
6.11 and Theorem 6.12]. 

Corollary 5.3 (combinatorial product). Introduce a product on the sub space CP^ 
by setting 

(5.8) \*fi:=Sx{A)fi, \/\fieP+. 
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Then forz = l (CP+, *) is a unital, associative and commutativ^ algebra isomorphic 
to the fusion or Verlinde algebra -F^ j^- Moreover, the renormalised Bethe vectors h\ = 
b\/ {h\,h\)^ are idempotents with respect to this product, i.e. 

(5.9) bx*b^ = 5x,J)\, yX, n eV<n-i,k ■ 

Thus, the completeness of the Bethe ansatz is equivalent to the semi- simplicity of the 
fusion algebra. 



Proof. The proof of the result ( 15 .Sp can be given in one hne, hence we repeat it 
here for the reader's convenience, 

/ ,\ \ba:l^) I Sqf(jS u,(jS\(j , ST^ I S\(jS uLcfSua \ 

'^(-^)^ = E M = E = E (E ) " • 

We aheady pointed out above that ( 15 .Sp imphes that (60-, A) = Sxa/S^^ and {b^^hcr) = 
\S%a\~'^ ^ whence b^j = S^j^cr Tlx It now follows that 

S 

bp*ba = S(lip S\pS\{A)b^ = — ^ SxpSxaba = 5a, pK , 
A "^0- A 

where in the last step we have used unitarity, 5" ■ 5"* = 1, of the modular S- matrix. □ 



§ 5.1. Fusion matrices as afRne plactic Schur polynomials 

It is well-known that the fusion matrices N^'^^ := (Af-'^^'^) - p+ form a representa- 

A A/x 'J^2!l k 

tion of the fusion ring. From the existence of the eigenbasis ( 15. 3 p and ( 15. 8p one deduces 
the next corollary which states that the affine plactic Schur polynomials ( 15. 7p . when 
restricted to the subspace CP^, are identical with the fusion matrices. 

Corollary 5.4. Denote by T'^ ^ C End(CP^) the subalgebra generated by the 
(restricted) affine plactic Schur polynomials {sx{a)k}xev<n-i k- ''^^P ^x{'^)k ^ [sa*] 
provides an isomorphism (z = 1), T'^ = C[ei, . . . , ek]/T(¥n,k) ■ In particular, one has 

(5.10) sx{a)kSf,{a)k = J^l^^''' s^{a)k ■ 

i^e'P<„-i,k 

Remark. It is common knowledge within the statistical mechanics community 
that a set of commuting transfer matrices is the distinguishing property of an integrable 

•^For arbitrary z the product is still associative but ceases to be commutative. This is different from 
(151 Theorem 6.12, eqn (6.33)] where the product was defined in terms of s^{A) with A being the 
partition obtained from A by adding k — Xi columns of height n. This introduces an additional 
^-dependence which renders the product commutative. 
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or exactly solvable lattice model [T]. Due to the development of the quantum inverse 
scattering method by the Faddeev school, it is the noncommutative structures, the Yang- 
Baxter algebras discussed in Section 3, which have been the centre of attention. The 
result ( 15.1 op shows that also the commutative (sub)algebra, the "integrals of motion" , 
have an interesting structure with applications in representation and, here, conformal 
field theory. 



§ 6. Recursion formulae for fusion coefficients 

The quantum mechanical interpretation via the Hamiltonian (13.241) identifies the 
fusion ring as the fc-particle superselection sector of the state space CP"*". The physical 
picture of creating and destroying particles via the maps (13. 3p and (13. 4p suggests to 
investigate how fusion coefficients at different levels are related. 

Let us start from the simple observation that according to (I3.35P the operator A{u) 
does not depend on an- We thus find from (13. 7p that A{u)lPi = (piA{u) and, hence, 

(6.1) if^T{u)^l = + z^^ip^A{u)^l = T{u) . 

Similarly, we can argue that Qo,o(^)'/'n = ^nQo,o{'^) (I3.39P and exploiting once more 
(13. 7p we arrive with the help of (I3.13P at 

(6.2) ^nQin)^l = Q{n) . 

We stay in the particle picture and set ^ = 1. Then the action of hk{A)\^p+ is par- 
ticularly simple: all particles on the su(n) Dynkin diagram are shifted by one position, 
hk[A)m = (m^i, mi, m2, . . . , mTj,_i). Because of (I3.2ip . the above commutation rela- 
tions (16. ip and (16. 2p generalise for 2; = 1 to the maps with 1 < z < n. If we recall 
that T and Q are the generating functions of the affine plactic elementary and complete 
symmetric polynomials we obtain immediately the following: 



Proposition 6.1 (recursion formulae). For any /i, z> G we have the identi- 



ties 



where < r < k and i = 1, 2, . . . , n. 



Example 6.2. Set n = 5 and fj, = (2, 2, 1, 0). Then for r = 3 we find with help 
of the algorithm of Section 4.1 the following expansion at level k = 2, 
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To verify the first identity in ( 16 .Sp for i = 1 observe that the first two terms in the 
foUowing expansion for = (3, 2, 1) at level /c = 3, 



+ 



+ 



+ 



are obtained from the k = 2 expansion by adding a one-column. 



Another set of relations follows from the recursion relations (16.11) and (|6.2p . As we 
discussed earlier, the boundary parameter z allows us to project on the finite plactic 
algebra A' by setting z — 0. With the help of (13.360 . (I3.4ip and (15. 8p one proves the 
following statement. 

Proposition 6.3. Let /i, z> G and /x, u the corresponding partitions under 



( li.^p . then we have for r = 0, l,...,7i — 1 the following expression in terms of Littlewood- 
Richardson numbers 



(6.4) 



(i'-)A 



1, 



if either = ui and u/ ^ = (1' 



or + 1 = ui and (fi*^i'/(fi\iJi = {V ^) 
0, else 



In case of an horizontal strip of length r = 0,1,..., /c one has instead the recursion 
relation 



(6.5) 



^,(k)u I 1' if ixi = vi andv/ix= {r) 
■^(7)A = 1 A^^^zi)'^^^ else 



Example 6.4. Again we set n = 5 and (i = (2, 2, 1, 0). Then at level k = A one 
computes via the above algorithm the expansion 



+ 



+ 



Since for all partitions v appearing in the expansion one has ui > only the second 
case in (16. 5p applies and, thus, we find the following nonzero fusion coefficients at level 
/c = 3, 

.r(3), (2,2,2,1) _ «r(3), (3,2,1,1) _ .r(3), (3,2,2) _ . 
■^^(2) (2,2,1) -■'^(2) (2,2,1) --^^(2) (2,2,1) " ^ ' 

The latter can be verified by using again the presentation ( 14. 8p in the ring of symmetric 
functions and the resulting algorithm or the Verlinde formula. 
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Remark. Setting z = in ( 14 ■4p it appears at first sight that we should expect to 
obtain the cohomology ring of the Grassmannian, 

H* {Grn-i,n+k-i) — ^[ei, . . . , e^]/ (/in, ■ ■ • , hn+k-l) ■ 

However, because of the ^-dependence entering the construction of the eigenvectors 
(15 -Sp one cannot conclude that the product (15 .Sp specialises to the cup product in 
H* {Grn-i,n+k-i) ■ Hence, the structure constants of the algebra defined via (15.80 be- 
come in this special case Littlewood-Richardson coefficients with the additional con- 
straint that = vi. 

§ 7. Conclusions 

While we have focussed here on a particularly simple integrable model and a very 
special ring, the observations made should generalise to a wider class of exactly solvable 
lattice models. Let us outline the general concept. 

Starting point in the construction of exactly solvable lattice models are in gen- 
eral solutions to the Yang-Baxter, or more generally, the star-triangle equation. The 
overwhelming majority of these solutions can be constructed with the help of repre- 
sentations of some noncommutative algebras, such as g-deformed enveloping algebras 
of Kac-Moody algebras (Drinfel'd-Jimbo "quantum groups") or elliptic generalisations 
thereof. For the example at hand the noncommutative algebras in question are the 
phase and affine plactic algebra and it has been explained in [16] how these have their 
origin in the g-deformed enveloping algebra UqSl{2). Once the solutions are known ex- 
plicitly one can interpret their matrix elements as the Boltzmann weights of a statistical 
vertex model - depending on a free parameter - and construct the corresponding trans- 
fer matrices to compute its partition function. Because the Boltzmann weights satisfy 
the Yang-Baxter equation the transfer matrices commute, thus defining a commutative 
(and associative) algebra or ring despite being built from the generators of a noncom- 
mutative algebra. We have seen how the transfer matrices are generating functions for 
polynomials in the alphabet (ai, . . . , a^) and while the letters of this alphabet do not 
commute the affine plactic Schur polynomials (15.70 do. 

Via the Bethe ansatz one then computes the idempotents of this commutative 
algebra, showing that it is semi-simple, and also its structure constants, the fusion 
coefficients and their expression in terms of the Verlinde formula. While the algebraic 
Bethe ansatz employed here is special to the t/qSl(2) case, generalisations of it which 
are applicable to higher rank, such as the nested, analytic or coordinate Bethe ansatz 
might be used instead. It is true for the majority of models solvable by the Bethe 
ansatz that the Bethe ansatz equations determining the eigenvectors of the transfer 
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matrix (or the idempotents of the commutative algebra) are given in terms of some 
elements in a commutative polynomial ring TZ and hence define an affine variety V. In 
the example discussed here the polynomial ring TZ has been the the ring of symmetric 
functions TZ — C[ei, . . . , e/s]. In the next step one determines the nuUstellen ideal I(V) 
and considers the quotient TZ/Xiy) which is isomorphic to the ring generated by the 
transfer matrices, namely we saw in the model discussed here that the eigenvalues of the 
affine plactic Schur polynomials for fixed particle number or level k could be identified 
with elements in the quotient ring; see Corollary 15.41 

In general it is not true that the Bethe ansatz equations can be solved, i.e. the 
variety V is not known explicitly. Nevertheless one might still be able to determine the 
nullstellen ideal X(V) and perform computations in TZ/I(V), similar to the computations 
of the fusion coefficients performed in Section 4, where we only used the abstract from 
of the polynomial equations (14. 4p but not the explicit solutions ( 15. ip . 

From these observations a natural classification question arises: can the commuta- 
tive algebras arising from integrable vertex models associated with the quantum groups 
UqQ be identified and do their structure constants have a similar representation theoretic 
interpretation? We hope to address this question in future work. 

Acknowledgments. The research of the author is carried out under a University 
Research Fellowship of the Royal Society. Part of the results summarised where pre- 
sented during the workshop Infinite Analysis 10 ''Developments in Quantum Integrable 
Systems''^ held at the Research Institute for Mathematical Sciences, Kyoto, Japan during 
June 14-16, 2010. The author wishes to thank the organisers. Professors Michio Jimbo, 
Atsuo Kuniba, Tetsuji Miwa, Tomoki Nakanishi, Masato Okado, Yoshihiro Takeyama 
for their kind invitation and hospitality. 

References 

[1] Baxter R. J., Exactly solved models in statistical mectianics, Academic Press Inc, Harcourt 

Brace Jovanovich Publishers, London, 1989 (Reprint of the 1982 original). 
[2] Bogoliubov N. M., Izergin A. G., and Kitanine N. A., Correlation functions for a strongly 

correlated boson system, Nuclear Phys. B, 516 (1998), 501-528. 
[3] Bogoliubov N. M., Izergin A. G., and Korepin V. E., Quantum inverse scattering method 

and correlation functions, Cambridge Monographs on Mathematical Physics, Cambridge 

University Press, Cambridge, 1993. 
[4] Di Francesco P., Mathieu P., and Senechal D., Conformal field theory, Graduate Texts in 

Contemporary Physics, Springer- Verlag, 1997. 
[5] Fomin S. and Greene C, Noncommutative Schur functions and their applications. Discrete 

Math., 193 (1998), 179-200. Selected papers in honor of Adriano Garsia (Taormina, 1994). 
[6] Fulton W., Young tableaux, volume 35 of London Mathematical Society Student Texts, 

Cambridge University Press, 1997. 



The su(n) WZNW fusion ring as integrable model 



33 



[7] Gepner D., Fusion rings and geometry. Comm. Math. Phys., 141 (1991), 381-411. 
[8] Gessel I.M. and Krattenthaler C, Cylindric Partitions, Trans. Am. Math. Soc. 349 (1997), 
429-479. 

[9] Goodman F.M. and Nakanishi T., Fusion algebras in integrable systems in two dimensions, 

Phys. Lett. B, 262 (1991), 259-264. 
[10] Goodman F.M. and Wenzl H., Littlewood- Richardson coefficients for Hecke algebras at 

roots of unity. Adv. Math., 82 (1990), 244-265. 
[11] Krattenthaler C., Guttmann A. and Viennot X., Vicious walkers, friendly walkers and 

Young tableaux: II. With a wall Jowrna/ of Physics A, 33 (2000), 8835-66. 
[12] Krattenthaler C., Guttmann A. and Viennot X., Vicious walkers, friendly walkers, and 

Young tableaux. III. Between two walls. (English summary) Special issue in honor of 

Michael E. Fisher's 70th birthday (Piscataway, NJ, 2001), J. Stat. Phys., 110 (2003), 

1069-1086. 

[13] Kac V. G., Infinite- dimensional Lie algebras. Cambridge University Press, second edition, 
1985. 

[14] Kac V. G. and Peterson D. H., Infinite-dimensional Lie algebras, theta functions and 

modular forms. Adv. in Math., 53 (1984), 125-264. 
[15] Korff C. and Stroppel C., The sl(n)fe-WZNW fusion ring: A combinatorial construction 

and a realisation as quotient of quantum cohomology. Adv. Math., 225 (2010), 200-268. 
[16] Korff C., Noncommutative Schur polynomials and the crystal limit of the C/qSt(2)-vertex 

model, J. of Phys. A: Math, and Theo., 43 (2010) 434021. 
[17] Lascoux A. and Schiitzenberger M.-P., Le monoi'de plaxique. In Noncommutative 

structures in algebra and geometric combinatorics (Naples, 1978), volume 109 of Quad. 

"Ricerca Sci." (1981), 129-156. 
[18] Macdonald I. G., Symmetric functions and Hall polynomials. Oxford Mathematical Mono- 
graphs, The Clarendon Press Oxford University Press, second edition, 1995. 
[19] Verlinde E., Fusion rules and modular transformations in 2D conformal field theory. 

Nuclear Phys. B, 300 (1988) 360-376. 
[20] Walton M. A., Fusion rules in Wess-Zumino-Witten models. Nuclear Phys. B, 340 (1990), 

777-790. 



